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We report on the latest results on the running coupling of two flavour QCD in the Schrodinger functional 
scheme. Results for the step scaling function are obtained from simulations on lattices L/a = 8 and L/a = 16 
which confirm the first results from lattices L/a = 4,5,6 presented one year ago by the ALPHA collaboration. We 
also discuss some algorithmic aspects, in particular concerning the occurrence of metastable states. A modified 
sampling, in order to estimate the proper weight of these states in the path integral, is proposed and tested. 



1. INTRODUCTION 

The relation between the short distance regime 
of QCD and its low energy sector is a problem 
which can be addressed by lattice techniques. 
The A parameter, characterizing the coupling at 
large energy, can be related to a hadronic quan- 
tity by a recursive finite-size technique, avoiding 
in such a way the multiple scale problem. This 
method has been successfully applied to the zero 
flavour approximation of QCD [|] and first results 
extending it to QCD with two massless flavours 
have been presented last year by our collabora- 
tion 1^ . Here we report about more recent results 
which allow for a better estimate of systematic ef- 
fects on our determination of the A parameter in 
the two flavour theory. This result is essential also 
in the computation of the running quark mass Q . 

We consider the Schrodinger functional (SF), 
defined as the Euclidean partition function of 
QCD on a cyhnder of size x T 



with periodic (up to a phase) boundary con- 
ditions in the spatial directions and Dirichlet 
boundary conditions at xq — 0, T . The spatial 
links at xq = Q,T are fixed to diagonal SU(3) 
matrices specified in terms of L and one angle 
rj while the quark fields are given by Grassmann 
values p, p and p' , p' , which are used as sources 
and then set to zero. 

We work with an 0(a) improved Wilson ac- 
tion and improved operators, employing the non- 
perturbative csw ID coefficient and one and two- 
loop values for the boundary coefficients Ct and 
Ct [|J^. The improvement coefficient ca is also 
set to its one-loop value. It enters our definition 
of quark mass m via the PCAC relation. The 
mass is required to be consistent with zero [G| . 

For T = i, the SF couphng g'^ is defined [m,0| 
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with k fixed by the condition = gi + 0((7q). 
Here L plays the role of a renormalization scale. 

2. NUMERICAL RESULTS 

The main quantity in the present computation 
is the step scaling function (SSF) 
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= 5^(2i)|g2(L)^„^„=o , (3) 

which can be viewed as a non-perturbative in- 
tegrated form of the /3-function. Once a{u) has 
been computed in the continuum hmit and for a 
large enough range of values, the equation 



a{fiL/2))^fiL) 
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can be recursively solved n times (for explicit 
equations see e.g. ^) starting from the hadronic 
scale imax with g^(Lmax) = "max- In this way 
we obtain values for g^(2^"Lijiax)- For large n 
this coupling is perturbative and we compute the 
A parameter in units of Lmax making use of the 
three-loop f3 function in the SF scheme 
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Later we will have to relate Lmax to a physical 
scale, e.g. by determining Lmax-pTr- Our new sim- 

L/a u E M E 

8 0.9807(17) 1.0745(55) 1.508(4) 1.716(14) 

8 1.1818(29) 1.3338(58) 2.014(10) 2.475(31) 

8 2.479(13) 3.348(48) 

4 3.334(11) 5.513(42) 

6 3.326(20) 5.62(9) 

Table 1 

New data for the lattice SSF Yj{u,a/L). ct((?o) 
was set to its one-loop value for the left part, 
and to two-loop for the two rightmost columns. 

ulations and results are summarized in Table 1. 
These allow for more reliable continuum limit ex- 
trapolations of the SSF. 

We extract the continuum limit values for the 
SSF using data for the two-loop improved observ- 
able S(2)(u,a/L) § from L/a = 6 and L/a = 8. 
The scaling of E^^^ vs. (a/L)"^ is shown in Fig- 
ure |l| for two different values of u. For compar- 
ison we included in the plot the continuum re- 
sults in the quenched case (black squares) for the 
same value of the couplings. The function a(u) 
is obtained interpolating the continuum limit es- 
timates by a sixth-order polynomial with the 
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Figure 1. Scaling plot for YP'^ . Black squares 
indicate the continuum limit extrapolation for 
TVf = 0. 

first three coefficients constrained by perturba- 
tion theory. Starting from Umax = 3.3 or Mmax = 5 
we can then estimate the quantity ALmax as dis- 
cussed above. Results are summarized in Table 2. 
We use the difference between continuum values 
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6/8 L/a = 5 
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1.82(5) 


1.87 
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1.23(5) 1.27 
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1.84(6) 


1.89 


7 


1.25(6) 1.28 


7 


1.85(7) 


1.91 


8 


1.26(7) 1.30 



Table 2 

Values for — ln(AL,nax) for two choices of Lmax- 

and data from L/ a — bio estimate our systematic 
uncertainty, which we linearly add to the statis- 
tical errors. This yields 



ln(AL„iax) 
ln(AL;,aJ 



4.85(13) 
-1.26(11) 



max ) 
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Comparing with the results in Q, we see that 
the systematic uncertainty on In(ALmax) is nearly 
halved, while the central values are fully consis- 
tent. The non-perturbative evolution of Oi{^) = 
(p{L)/A-K = 1/L) starting from cp = 5 is 
plotted in Figure ||. Statistical and systemat- 
ical errors (difference between continuum limit 
and L/a — 5) are too small to be visible in the 
plot. The figure shows that differences between 
Monte Carlo results and three-loop perturbation 
theory become appreciable for large values of the 
coupling, while for high energies the perturbative 
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approximation very closely reproduces the non- 
perturbative results. 
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Figure 2. Running of a — g'^/An in the SF scheme. 

3. METASTABLE STATES 

For the largest values of the coupling in our 
study {g^ ~ 3.3 and g"^ ~ 5.5) we observed a 
long tail towards negative values in the distribu- 
tion of dS/drj. In addition, measuring by a cool- 
ing procedure the pure gauge contribution to the 
action Sq°\ and to the coupling dS^"^ / dr], we 
observed occurrences of metastable states in our 
Monte Carlo histories. The action S(?°^ for these 
states is consistent with the value for a secondary 
solution of the field equations , given our choice 
for the boundary fields. This solution appears to 
be a local minimum by numerical evidence. 

In order to properly estimate the weight 
of these states we enhanced their occurrence 
through a modified sampling adding to the 
HMC effective action a term 
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where 70 and w-f are fixed values, while 7 is a dy- 
namical variable. The expectation values in the 
original ensemble are then obtained by reweight- 
ing properly the observables. 

Defining a quantity q whose value is 1 for 
metastable states and otherwise (so that {5qi) + 
(Sqo) = 1)) for an observable O one can write 



(O) = {SqlO) + (SgoO) 

= (V0) + (0)i(l-(<5,i)), 



(8) 



where {0)i = {SqoO) / {Sqo) . If the main contri- 
bution to (O) comes from the q — sector then 
a precise estimate of (O) just requires a precise 
estimate of (0)i, which can be obtained by an 
algorithm which samples only the q = sector, 
together with rough estimates of {5qi) and {SqiO) 
which can be obtained by the modified sampling. 

By a set of simulations on L/a — 8 and 12 
using this approach we could estimate the effect 
of metastable states on g"^ to be nearly 0.10(2)% 
for g^ ~ 3.3. This result is independent, within 
errors, from L/a. Moreover, for this value of 
the coupling, we noticed a slightly better behav- 
ior for the PHMC algorithm in comparison with 
HMC (smaller values of 70 needed in order to gen- 
erate a certain number of transitions from/to a 
metastable state). 

For g^ ~ 5.5, the effect of metastable states is 
much larger, as expected, resulting in a 4 to 5% 
contribution. However their sampling is also en- 
hanced in this case using both PHMC or HMC 
already in the original ensemble. Indeed we re- 
peated the L/a = 12 simulation for g^ ~ 5.5 using 
ordinary PHMC as in Ref. ||^ but disentangling 
the occurrence of metastable states. This turned 
out to be around 6% and the result for g^ was 
fully consistent with the number in Ref. 1^. 
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